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^ ! Abstract 

Binary optimal codes often contain optimal or near-optimal subcodes. In this 
'^ ■ paper we show that this is true for the family of self-dual codes. One approach is to 

r~^ ! compute the optimum distance profiles (ODPs) of linear codes, which was introduced 

CN ' by Luo, et. al. (2010). One of our main results is the development of general algorithms, 

called the Chain Algorithms, for finding ODPs of linear codes. Then we determine 
the ODPs for the Type II codes of lengths up to 24 and the extremal Type II codes 
of length 32, give a partial result of the ODP of the extended quadratic residue code 
CN ' Q48 of length 48. We also show that there does not exist a [48, k, 16] subcode of q^g 

for k > 17, and we find a first example of a doubly-even self-complementary [48, 16, 16] 
code. 
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1 Introduction 

One of the main problems that has arisen in Coding Theory is the search for optimal codes 
with the largest size given a minimum distance or optimal codes with the largest minimum 
distance given a size [T2| 121] ITf] . There has been extensive work in this direction [9j . Some 
well-known families of codes, such as the Reed-MuUer codes or the cyclic codes, contain 
notable subcodes. However, comparatively little attention has been paid to the subcodes of 
an optimal linear code in general. It is a natural concern to determine which linear codes 
contain optimal (or near-optimal) subcodes. Among linear codes, we suggest self-dual or 
self-orthogonal codes since their possible non-zero weights jump by 2 or 4. Thus there is a 
possibility to get subcodes with a large minimum distance. 

In fact, self-dual codes have been one of the most active topics in algebraic coding 
theory since V. Pless |20] started to classify binary self-dual codes in 1972. These codes 
have interesting connections to groups, t-designs, lattices, and theta series [121 [13 l2l]- 
Furthermore, many extremal self-dual codes often turn out to be the best among the linear 
codes with the same parameters. Nevertheless, little attention has been paid to the subcodes 
of self-dual codes. 

We plan to construct optimal (self-orthogonal) subcodes of a given linear (self-dual) code. 
In order to construct finite-state codes, PoUara, Cheung and McEliece [23] constructed the 
first [24, 5, 12] subcode of the binary Golay [24, 12, 8] code, improving a previously known 
[24, 5, 8] subcode. Maks and Simonis [18] have shown that there are exactly two inequivalent 
[32, 11, 12] codes in the binary Reed-Muller code R{2, 5) which contain -R(l, 5) and have the 
weight set {0,12,16,20,32}. 

We show that in many cases optimal subcodes can be obtained by computing optimum 
distance profiles (ODPs), a concept introduced by Luo, Han Vinck, and Chen [16]. The 
authors [16] considered how to construct and then exclude (or include, respectively) the basis 
codewords one by one while keeping a distance profile as large as possible in a dictionary order 
(or in an inverse dictionary order, respectively). Thus fault-tolerant capability is improved 
by selecting subcodes in communications and storage systems. The practical applications 
are found in WCDMA [TT], [26] and address retrieval on optical media [2]. 

In [5] and 116], the authors give results on the ODPs of the binary Hamming [7, 4, 3] code, 
the binary and ternary Golay codes, Reed-Solomon codes, the first-order and second order 
Reed-Muller codes. Recently, Yan, et. al. [2H] considered the optimum distance profiles 
of some quasi-cyclic codes and proposed two algorithms, called the "subcodes traversing 
algorithm" and "supercodes traversing algorithm." These algorithms enumerate all subcodes 
of a given code. Hence they are rather inefficient in finding ODPs of linear codes with a 
relatively large dimension. Their examples have dimension 10 only. Hence we ask the 
following two questions. 

(i) Is there an interesting class of linear codes whose ODPs are not known yet? 
(ii) Is there an efficient algorithm to compute ODPs of linear codes? 

For question (i), we choose a class of self-dual codes since the structure of these subcodes 
is surprisingly less known. For question (ii), we propose two full algorithms based on cosets, 
called the Chain Algorithms and two random algorithms to find ODPs of the codes. These 
algorithms look at a chain of subcodes of a given code and consider the equivalence of 



the codes with the same dimension. Hence they are more efficient than the subcodes and 
supercodes traversing algorithm 



In this paper, we give the ODPs of Type II self-dual codes of lengths up to 24 and the 
five extremal Type II codes of length 32, give a partial result of the ODP of the extended 
quadratic residue code ^43 of length 48. Moreover, we show that each of the five Type 
II [32, 16,8] codes contains the two optimal [32, 11, 12] codes, which was previously known 
only for the Reed-Muller code -R(2, 5). We also construct a [48, 14, 16] code and an optimal 
[48, 9, 20] code from the extended quadratic residue code q^s oi length 48. Both codes are 
not equivalent to the best known codes of the same parameters in the Magma database |1]. 
We also show that there does not exist a [48, k, 16] subcode C of ^43 for k > 17. We find 
a first example of a doubly-even self-complementary [48, 16, 16] code. Such a code was 
previously not known to exist. Only one singly-even self-complementary [48, 16, 16] code 
was found by A. Kohnert [12]. Similarly we construct [72,29,16], [72,23,20] codes which 
are not equivalent to the best known codes. Further we construct a new self-orthogonal 
[72,35, 16] code with Aiq = 129972. All our computations were done using Magma [1]. 

2 Preliminaries 



We refer to p2] for basic definitions and results related to self-dual codes. All codes in this 
paper are binary. A linear [n, k, d] code C of length ra is a fc-dimensional subspace of GF{2)^ 
with the minimum (Hamming) weight d. Two codes over GF{2) are said to be equivalent 
if they differ only by a permutation of the coordinates. The dual of C, denoted by C"*" is 
the set of vectors orthogonal to every codeword of C under the Euclidean inner product. If 
C = C"*", C is called self-dual. A self-dual code is called Type II (or doubly-even) if every 
codeword has weight divisible by 4, and Type I (or singly-even) if there exists a codeword 
whose weight is congruent to 2 (mod 4). 

Let C be a binary self-dual code of length n and minimum distance d{C). Then d{C) 
satisfies the following 



diC) < <! 4 



n 

.24 
n 

24 



+ 4, if n 7^ 22 (mod 24), 
+ 6, if n = 22 (mod 24). 

A self-dual code meeting one of the above bounds is called extremal. 

A subcode C of a linear code C with minimum distance d' = d{C') > d{C) is maximal 
if there is no subcode C" such that C C C" C C and d{C") = d'. Given d' > d{C), the 
maximum of the dimensions of maximal subcodes C with d{C') = d' is called the maximum 
dimension with respect to d' . Given n and fc, a hnear [n, fc, d] code is minimum distance 
optimal if d is the largest possible. (Grassl's online table [9] is a good source for reasonable 
lengths and dimensions for finite fields of order up to 9.) Given n and d, a linear [n, fc, d] 
code is dimension optimal if k is the largest possible [121 P- 53]. We raise the following 
natural question. Given a binary self-dual code C and any non-zero weight d' > d{C), how 
do we find a subcode with maximum dimension with respect to d'7 In general, this question 
seems very difficult since theoretically we should know all subcodes. On the other hand, 
there has been another approach related to this problem, as described below. 



Let C be a binary [n, k] code and let Cq = C . A sequence of linear subcodes of C, 
Cq ^ Ci ^ ■ ■ ■ ^ Ck-i is called a subcode chain, where the dimension oi Ci is k — i for 
i = 0, ... A; - 1. (If we let Ck = {0} and Vi = Ck-^{i = 0, ■ ■ ■ , A;), then {0} = K) C Vl C 
• • ■ C Vfc = C is known as a complete flag [19j.) 

Let (ij = (i(Ci) be the minimum distance of Ci. Then the sequence do < di < ■ ■ ■ < dk-i 
is called a distance profile of C (see [5], [T6] for details). A generator matrix such that its 
first k — i rows (i.e., the remaining rows after removing its i rows from the bottom) form 
a generator matrix of Cj for < i < A; — 1, is called a generator matrix with respect to the 
distance profile. 

For any two integer sequences of length k, a = uq, . . . , ak-i and b = bo, ... , bk-i, a is 
called an upper bound on b in the dictionary order if a is equal to b or there is an integer t 
such that 

Oj = bi for < z < t — 1, and at > bf. 

On the other hand, a is called an upper bound on b in the inverse dictionary order if a 
is equal to b or there is an integer t such that 

Qi = bi for t + l<i<k — 1, and at > bf. 

A distance profile of the linear block code is called the optimum distance profile (or ODP 
for short) in the dictionary order, which is denoted by ODP'^*'^[C](0), 0DP'^*'^[C](1), . . . , 
QDpd.^'^^C]{k — 1) if it is an upper bound on any distance profile of C in the dictionary 
order. Similarly, a distance profile of the linear block code is called the optimum distance 
profile (or ODPioi short) in the inverse dictionary order, which is denoted by ODP*"^[C](0), 
0DP*™[C](1), . . . , ODP*"^[C](/i; — 1) if it is an upper bound on any distance profile of C in 
the inverse dictionary order. To simplify notations, for a given [n, k] code C we may use 
Oi3prfic[q_ OI)P'^'^[C]{k-i) (resp. ODP*™[C]i= ODF''''"[C]{k-i)) so that we may easily 
interpret the corresponding subcode parameters: [n,i, ODP'^*'^[C]j] (resp. [n,i, ODP*™[C]i]). 
We also use ODP[C] to denote the optimum minimum distance profile in both orders. Note 
that for a given [n, k] code C over GF{q), the number of its subcode chains [16J is 

ngiM-ii=ri(^. 

t=2 t=2 ^ 

where Q[t,r] is the g-ary Gaussian binomial coefficient n^o ''^-jZi - Hence for large dimen- 
sions it will be very difficult to determine ODP of a linear code by a brute-force search. 

3 Relation between ODP and the maximum dimension 

The ODP of a code and the maximum dimension with respect to a minimum distance 
are related concepts. Note that the first minimum distance d' to appear in the ODP in 
dictionary order corresponds to a maximal subcode with maximum dimension corresponding 
to d'. However, after this term, maximal subcodes in the subcode chain do not necessarily 
imply the maximum dimension. This is an observation which follows from the definition of a 
maximal subcode and the definition of ODP; we formalize the theory in the following results. 



However, note that given a dimension k' < k there may be multiple minimum distances d' 
with respect to which k' is the maximum dimension. Therefore for the first proposition we 
define dk' to be the maximum of such minimum distances. 

Proposition 3.1. Let C he an [n, k] code. Let k' < k be given. Define dy = ina.x{{d' : k' is 
the maximum dimension in C with respect to d'}) and define dgpt to be the optimal minimum 
distance attained among all [n, k'] codes (many values available at JMl), then 



dic\ 



Proof. The claim dopt > dk' is clear since dopt is the maximum minimum distance possible 
among all [n, k'] codes. By the definition of dk', if C contains an [n, k', d'] subcode, then dk' > 
d'. Since ODP'^*^[C]fc. (respectively ODP*™[C]fcJ corresponds to a dimension ki subcode in 
the subcode chain having minimum distance ODP'^*^[C]fc^ (respectively ODP*"^[C]fc.), the 
preceding claim proves the proposition. D 

Corollary 3.2. Let C be an [n,k] code. Let k' < k be given. Define dy and dopt os above. 
If ODF'^-[C]k' = dopt or ODP--[C]k' 



dopt, then 



dici 



dopt = dk' = m^x{{ODP^''^[C]k', ODP'"'[C]k'}). 



The necessity of defining dk', in Proposition 13. ![ as a maximum is due to the fact that 
there may be multiple minimum distances yielding the same maximum dimension. An 
example where this occurs is the following: 



Example 3.3. Let C be the [6,3,1] code with the following generator matrix: 



G 



The maximum dimension with respect to (ii = 4 is 2, due to the fact that the first two rows 
of G generate a [6,2,4] subcode of C with the following generator matrix: 



1 11 


00 


1 00 


11 


00 


00 



Gi 



11 11 00 
11 00 11 



Similarly, the maximum dimension with respect to (i2 = 3 is 2; this is obtained by adding 
the third row of G to each row in Gi which yields a [6, 2, 3] subcode of C with the following 
generator matrix: 

" 01 11 00 " 
01 00 11 



G, 



Notice that in Proposition l3.1l we fix the dimension k'; a dual statement where we instead 
fix the minimum distance is the following. 



Proposition 3.4. Let C he an [n, k] code and let < j < k — 1. Suppose dj is a minimum 
distance appearing as ODP^^^[C]j or ODP™[C]j. Define kj to be the maximum dimension 
with respect to dj, then kj > j . 

Proof. The proof follows directly from the definition of maximal dimension with respect to 
dj, since a subcode with this maximal dimension will have dimension kj which is an upper 
bound on the dimension of any [n, j, dj] subcode. D 

The following proposition is a special case of Proposition 13. 4j this proposition states that 
in fact the first minimum distance in the dictionary order ODP corresponds to a maximal 
subcode with respect to that minimum distance. 

Proposition 3.5. Let C be an [n,k,d] code. Suppose that for some j , ODP^'^^[C]j is the first 
term in ODP greater than d. Then j is the maximum dimension with respect to ODP^^'^[C]j. 

Proof. If ODP'^*'^[C]j is the first term in ODP greater than d, then ODP'^*^[C]j+i = d where 
< j < k. Suppose to the contrary that j is greater than the maximum dimension with 
respect to ODP'^*^[C]j, then there must exist an [n,j + 1] subcode with minimum distance 
ODpdic^Q^, r|nj^ig implies ODF'^'^[C]j+i = ODP'^*^[C]j by definition of the dictionary order. 
Compiling this information we obtain the contradiction: d = ODP'^*^[C]j+i = ODP'^*'^[C]j > 
d. D 

Propositions I3.H I3.4[ and 13.51 give insight into the relation between maximum dimension 
subcodes and optimum distance profiles. If a code contains an optimal subcode (minimum 
distance optimal, dimension optimal, or both) there are many cases where this subcode 
appears in the subcode chain involved in an optimum distance profile. However, this is not 
always the case as in the following example: 

Example 3.6. Let C be the [6,5,1] code with the following generator matrix: 



G 



11 11 00 

11 00 11 

10 10 10 

10 10 00 

10 00 00 



By expurgating weight 1 vectors from C we may obtain [6, 4, 2] subcodes of C. Since 
there does not exist a [6,4,3] code (see [9]), we may conclude that 0DP'^*'^[C]4 = 2. By 
examining all [6, 4, 2] subcodes of C it can be determined that none contain a [6, 3, 3] subcode, 
and since no [6,3,4] code exists we obtain 0DP'^*'^[C]3 = 2. Finally, there is a unique 
[6, 2, 4] code (which has a single non-zero weight of 4); as this code is a subcode of at least 
one [6,4,2] subcode of C, and since there does not exist a [6,2,5] code we may conclude 
0DP'^*'^[C]2 = 4 and ODP*'^[C]i = 4. Therefore the optimum distance profile in dictionary 
order is ODP^^=[C] = [1,2,2,4,4]. 

Using similar arguments the ODP in inverse dictionary order is obtained as ODP*™ [C] = 
[1,2,2,3,5]. Notice that the first three rows of G generate an optimal [6,3,3] code (both 
minimum distance optimal and dimension optimal). Therefore the maximum dimension 

6 



with respect to ininiinum distance d' = 3 is k' = 3. However, the subcodes of dimension 
3 appearing in both ODP orders have minimum distance 2. An explanation for this phe- 
nomenon is that all supercodes of the [6,3,3] code in C have minimum distance 1. This is 
an example where equality is not possible in Proposition 13.11 and in Proposition 13. 4[ 



4 Algorithms based on cosets 

Given an [n, k, d] code C which has small length and dimension it may be relatively easy to 
examine its subcode structure by a brute force generation of all possible subcodes. However, 
as length and dimension increase this method becomes very time consuming; this is why we 
propose four algorithms based on cosets which are relatively efficient in comparison to the 
brute force search. The first two algorithms, called the Chain Algorithms are useful in the 
classification in the sense that when applying them we obtain a complete list of inequiv- 
alent subcodes (respectively supercodes), with prescribed minimum distance, contained in 
(respectively containing) the given code C; in this way, the redundant cases considered in 
a brute force search are eliminated. The remaining two Random Algorithms are random 
versions of the Chain Algorithms, and especially useful for very large length and dimension, 
where the exhaustive search is infeasible. The Random Algorithms can also give results 
much faster than the Chain Algorithms since all cases are not considered. 

(Subcodes) Chain Algorithm I: An algorithm to produce all maximal subcodes with 
maximum dimension k' and minimum distance d' > d. 

(i) Input: Begin with a binary [n, k, d] code D and a positive integer d' > d (such that 
there exists a codeword of weight d' in D). 

(ii) Output: Produce the maximum dimension k' among all maximal subcodes with min- 
imum distance d' and a list of inequivalent maximal subcodes of this dimension and 
minimum distance d'. 

(a) Initialize the set Bi = {-D"*"}. Begin with i = 1. 

(b) Build a set Bj+i of all inequivalent supercodes of dimension 1 higher of C for all 
C G Bj. In order to do this we add coset representatives from F^/C to each code 
C in Bi. 

(c) Check if d{C^) = d' for any code C G B^+i. If "No" for any C G B^+i, then 
repeat step (ii) by increasing i to i + 1. If "Yes", then output the maximum 
dimension k' = k — i + 1 and the set oi [n, k — i + 1, d'] subcodes of D. 

(Supercodes) Chain Algorithm II: An algorithm to find all [n, k, d] supercodes contain- 
ing an [n, k', d'] code with d' > d and k > k' 



(i) Input: Begin with a set Ck\d' of inequivalent [n, k', d'] codes (respectively self-orthogonal 
codes) with k > k' and d' > d. 



(ii) Output: For each code C in Ck',d'^ produce all [n, fc, d] codes (respectively self-orthogonal 
codes) containing C. 

(a) Begin by building a set of all inequivalent supercodes (respectively self-orthogonal 
supercodes) of dimension 1 higher of each code C in Ck'4' with minimum distance 
greater than or equal to d. In order to do this we add coset representatives from 
Fg/C (respectively C^ jC if C is self-orthogonal) to each code C in Qw.d' and 
keep a set of inequivalent supercodes C^'+i generated in this way. 

(b) Repeat the first step, by replacing 0^',^' with Cfc/+i until the set of inequivalent 
codes which are generated have dimension k. 

(c) Stop once dimension k is reached. For each code C in Q^'^' output all [n, k, d] 
supercodes of C . 

Analysis and comparison of our algorithms: 

Given an [n, /c] code C, the search for subcodes of dimension k' may be conceptualized 
as a search tree with root C and each node of branch distance h from C given by a [n, A; — 6] 
subcode. A brute-force search of the subcodes of dimension k' for an [n, A;] code searches 
through all branches of the search tree up to distance k — k'; this search has complexity 

k 



given by the Gaussian binomial coefficient 



k' 



The Chain Algorithms greatly reduce 



this search by "pruning" the search tree in two manners. First, we keep only inequivalent 

subcodes (resp. supercodes) at each branch level (in addition this keeps the search efficient 

memory- wise) . Second, branches can only extend from subcodes that were preserved in 

the previous step creating a chain of subcodes. In comparison, the algorithms given in 

Yan, et. al. |28j construct all subcodes of the same dimension not necessarily in chains of 

codes; this method corresponds to searching all nodes at a given branch distance (many of 

which are redundant). 

For example, a brute-force search of the subcodes of dimension k' for an [n, k] code 

r j^ ] 
has complexity given by the Gaussian binomial coefficient , , . In Section 15.41 for some 

L ^ J2 

[32, 16, 8] codes we determine the maximum dimension subcode with respect to (i = 12 to 
have dimension 11. A brute- force subcode search (such as the subcodes traversing algorithm 

^^ = 120, 843, 139, 740, 969, 555 subcodes; this task 
2 



in [2E]) would have to enumerate 
is not feasible. 



11 



Example 4.1. As a more concrete example, we determine the ODPs for the four optimal 
[28, 7, 12] self-complementary codes classified in [S]. These codes are doubly-even with non- 
zero weights 12,16,28. We begin with a [28, 3, 16] constant weight code (meaning the only 
non-zero weight is 16). There is only one such code due to the fact that all non-zero 
codewords must intersect in exactly 8 positions; if the first two basis vectors are fixed, then 
there is only one possibility (up to coordinate permutation) for the third basis vector. By 
adding the all-one vector to the constant weight code we obtain a [28, 4, 12] code with the 



following generator matrix: 



G 



[28,4,16] 



nil 0000 0000 nil oooo iiii iiii 

0000 nil 0000 nil nil oooo nil 

0000 oooo nil nil nil nil oooo 

nil nil nil nil nil nil nil 



Applying (Supercodes) Chain Algorithm II to this generator matrix (and only keeping 
doubly-even supercodes) we obtain all four self-complementary [28, 7, 12] codes with the 
following generator matrices: 



G^[28,4,16] 

0100010001011010010010111001 
0010011101110111001111001100 
0001000100011110010101010011 



^[28,4,16] 

0100010001001101001110101010 
0010011101110111001111001100 
0001000100011110010101010011 



G^[28,4,16] 

0100010001001011010100110110 
0010011101110111001111001100 
0001000100011110010101010011 



^[28,4,16] 

0101001100111010000010011010 
0011000000110011001100110011 
0000011001011010010110101100 



Let C be any [28, 7, 12] self-complementary code. Since the [28, 3, 16] subcode is optimal, 
in light of Lemma l3.5l we determine 0DP'^"^[C]3 = 16. As a [28,3,16] subcode cannot contain 
the all-one vector, we determine the ODP in dictionary order: 

ODP°'^"[C] = [12, 12, 12, 12, 16, 16, 16]. 

The ODP in inverse order is clear since any supercode of the repetition code, containing 
a weight 16 vector, must also contain a weight 12 vector. Hence 

Qj3pm,;^q = [12, 12, 12, 12, 12, 12, 28]. 

We now introduce the random algorithms which are random versions of the above coset 
algorithms: 

Random (Subcodes) Algorithm I: An algorithm to search for maximal subcodes 

(i) Input: A linear code C with parameters [n, k, d] and d' > d where A^/ is non-zero, 
(ii) Output: A maximal subcode C of C with d'. 



(a) Take any codeword x from C such that wt(x) > d'. Let Ci = (x). 

(b) Choose any coset representative y of C/Ci. Let Ci := {y) +Ci. Repeat this until 
d{Ci) = d'. 



(c) Repeat (b) until there is no coset representative such that d{Ci] 
C := Ci. 



d'. Let 



The below algorithm is somewhat opposite to Random Algorithm I. 
Random (Supercede) Algorithm II: An algorithm to search for codes containing good 
codes 

(i) Input: A (best known) linear code Ci with parameters [n, k, d] and d' < d. 
(ii) Output: A code C containing Ci with d' and k' > k (if such a C exists). 

(a) Let C:=Ci^. 

(b) Choose any coset representative y of C/Ci. Let Ci := {y) + Ci. Repeat this until 
rf(Ci) = d'. 

(c) Repeat (b) until there is no coset representative such that d{Ci) = d' . Let 
C ■= Ci. 

Example 4.2. Using their traversing algorithms, the authors [28] have determined ODPs 
of a quasi-cyclic [48, 10, 20] code C48 by finding all fc- dimensional subcodes of C which is 
extensive work. Using the above Random Algorithms, we have also computed ODPs of C48 
in the dictionary and inverse dictionary orders in a minute as follows: 

ODP'^*^[C48] = [20, 20, 20, 20, 24, 24, 24, 24, 32, 32], 
Qj^pint, ^^^^1 ^ pg^ 20, 20, 20, 20, 20, 20, 24, 28, 36] . 



5 ODP of Type II self-dual codes 

In this section, we determine the ODP of binary Type II codes of lengths up to 24 and the 
extremal Type II codes of length 32. 



5.1 



n = 



For length n = 8, there is a unique binary Hamming [8,4,4] code eg. It has two non-zero 
weights 4 and 8. It is clear that there is a unique subcode (1) of eg with ^4 = 8. Hence 

0DP[e8] = [4,4,4,8]. 

One generator matrix with respect to the ODP in the dictionary order is 

G{es) 



' 11 


11 


11 


11 " 


00 


00 


11 


11 


00 


11 


00 


11 


01 


01 


01 


01 



10 



5.2 



n 



16 



Next let us consider n = 16. There are two Type II [16,8,4] codes, denoted by rfig and 
2e8 [6] (blank represents 0): 



G{d 



Id) 



' 11 


11 








11 


11 








11 




11 






11 




11 






11 




11 






11 






11 




11 








11 


1 


1 1 


1 1 1 


1 


1 



and ^(268) = ^(es) © ^(es) 



The next higher weight in diQ is 8. We have constructed a [16, 5, 8] subcode of di^. This 
subcode is equivalent to the first order Reed-MuUer code i?(l,4) and hence is unique up 
to equivalence [27]. As there is no [16,6,8] code [9], we know that k = 5 is the maximum 
dimension with respect to d = 8. Since i?(l,4) contains the all-one vector, we have 

0DP[di6] = [4,4,4,8,8,8,8,16]. 

Considering some linear combinations of the rows of G{diQ), we give below one generator 
matrix with respect to the ODP in the dictionary order. 



G'id.e) 



' 11 


11 


11 


11 


11 


11 


11 11 " 


11 


11 


11 


11 








11 


11 






11 


11 




11 




11 




11 




11 
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1 1 


11 






11 








11 










11 




11 












11 _ 



In a similar manner, we have verified that 2e8 has a maximal [16, 5, 8] subcode, which is 
generated by the first five rows of G'{diQ). Hence we have 

ODP[2e8] = [4,4,4,8,8,8,8,16]. 
We give below one generator matrix with respect to the ODP in the dictionary order. 



^(268) 



■ 11 


11 


11 


11 


11 


11 


11 


11 


11 


11 


11 


11 










11 


11 






11 


11 






11 




11 




11 




11 




1 


1 


1 


1 


1 


1 


1 


1 










11 


11 














11 




11 












1 


1 


1 


1 
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As a summary, we have 
Theorem 5.1. 



0DP[di6] = ODP[2e8] = [4,4,4,8,8,8,8,16]. 



5.3 n = 24 



Consider n = 24. There are exactly nine Type II self-dual codes of length 24. These are de- 
noted by A2A{2di2), B2A{dw + 267), C2A{3ds), D24(44), ^24(^24), F24(6t/4), G2A{g24), die + 
eg, and Seg in the notations of [B], [22]. The first seven codes are indecomposable and the 
rest are decomposable. Note that G24((724) represents the binary Golay [24, 12, 8] code. 



PoUara, et. al. [23] constructed the first [24, 5, 12] subcode C24 of (724, improving a 



previously known [24,5,8] subcode. Note that C24 is unique [27], has only two non-zero 
weights 12 and 16, and has a [24, 2, 16] subcode C24 • As €24 satisfies the Griesmer bound, 
it has a generator matrix of which each row has weight 16 [27], [T2]. Hence it is easy to see 
that C24 is unique. 

Using this information, Luo, et. al. [16] have determined 

ODP*^[(724] ^ 
ODP*™[(724] ^ 



[8,8,8,8,8,8,8,12,12,12,16,16] 
[8,8,8,8,8,8,8,8,12,12,12,24]. 



However, less is known of the subcodes of the other Type II self-dual codes of length 
24. We have checked that the unique [24, 5, 12] code is contained in any of the nine Type II 
codes of length 24. 

Using (Subcodes) Chain Algorithm I we obtain inequivalent maximal [24, k', 8] subcodes 
of each Type II code of length 24 (with minimum distance 4). Then applying (Supercodes) 
Chain Algorithm II to the unique [24, 5, 12] code for each Type II code of length 24 (with 
minimum distance 4) we obtain a [24, k', 8] code equivalent to one of the maximal subcodes. 
Therefore we determine the ODP in the dictionary order of the Type II [24, 12, 4] codes as 
follows. 



5.2. 










ODI^''[2di2 


= 4,4,4,8,8,8,8,12,12,12,16,16 


ODP^'idio + 267 


= 4,4,4,8 


8 


8 


8,12,12,12,16,16 


ODP^'isds] 


= 4,4,8,8 


8 


8 


8,12,12,12,16,16 


ODF^'iAde 


= 4,4,8,8 


8 


8 


8,12,12,12,16,16 


ODI^'id24] 


= 4,4,4,4 


8 


8 


8,12,12,12,16,16 


ODP^'ied^ 


= 4,8,8,8 


8 


8 


8,12,12,12,16,16 


ODP'^'idie + es] 


= 4,4,4,8 


8 


8 


8,12,12,12,16,16 


ODP^'^lSes 


= 4,4,4,8 


8 


8 


8,12,12,12,16,16 



For each Type II [24, 12, 4] code we apply (Subcodes) Chain Algorithm I to the maximal 
[24, k', 8] subcodes (containing the all one vector) to obtain a [24, 4, 12] subcode (containing 
the all one vector). Therefore we may determine the ODP in the inverse dictionary order of 
the Type II [24, 12, 4] codes as follows. 
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Theorem 5.3. 



ODF''''[2di2] = [4,4,4 

OZ)/^™[dio + 2er] = [4,4,4 

OL'i^^isrfg] = [4,4,8 

ODF''''[Ade] = [4,4,8 

ODF''''[d24\ = [4,4,4 

0DF''''[Qd4 = [4,8,8 

ODF^''[d,e + es] = [4,4,4 

OZ)/^™[3e8] = [4,4,4 



8,12,12,12,24] 
8,12,12,12,24] 
8,12,12,12,24] 
8,12,12,12,24] 
8,12,12,12,24] 
8,12,12,12,24] 
8,12,12,12,24] 
8,12,12,12,24] 



Table [T] gives the maximum dimension with respect to minimum distance d for the Type 
II length 24 codes. 

Corollary 5.4. For each Type II length 24 code, there are maximum dimension subcodes 
with respect to d = 8, 12, 16, 24 (except 20 j that are involved in the subcode chain for 
the ODP in dictionary order or the inverse order. Furthermore, each Type II length 24 
code contains dimension optimal (and minimum distance optimal) subcodes with parameters 

[24, 5, 12], [24, 2, 16], [24, 1,24]. 



Table 1: Subcodes of All Type II codes of n = 24 



Codes 


max. dim. 


max. dim. 




with d = 8 


with d=12 


2di2 


9 


5 


dio + 2er 


9 


5 


3^8 


10 


5 


44 


10 


5 


d24 


8 


5 


6(^4 


11 


5 


di6 + eg 


9 


5 


Seg 


9 


5 


fl'24 


12 


5 



5.4 



n 



32 



As there are 85 Type II self-dual codes of length 32, we focus on extremal Type II self- 
dual [32, 16, 8] codes. There are exactly five Type II self-dual [32, 16, 8] codes, denoted by 
C81 (or 532), C82 (or rga, i?(2,5)), C83 (or 2^i6), C84 (or 8/4), C85 (I6/2) in the notation 
of [6] , [7] . Using symplectic geometric approach, Maks and Simonis |18] show that the second 
order Reed-Muller code r^2 contains exactly two inequivalent [32, 11, 12] codes, each of which 
further contains the first order Reed-Muller [32, 6, 16] code i?(l, 5). Note that any [32, 6, 16] 
code is equivalent to i?(l,5). Furthermore, Jaffe [13] proved using his language Split 
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that there exist exactly two [32, 11, 12] codes. These subcodes have optimal dimensions for 
each minimum distance. Hence Chen and Han Vinck [S] have determined the ODP in the 
dictionary order for r^2 as follows: 

ODP[r32] = [8, 8, 8, 8, 8, 12, 12, 12, 12, 12, 16, 16, 16, 16, 16, 32]. 

On the other hand, little was known of the subcodes of the other four extremal Type 
II [32, 16, 8] codes. We show that they also have the same optimum distance profiles as r32 
does. 

Using (Supercodes) Chain Algorithm II with Ck\d' = {-^(1)5)}, we independently con- 
struct two inequivalent [32,11,12] codes in r32 containing /2(1,5), denoted by RCi and 
RC2- We note that dim(i?Ci r\RC2) = 10. Using (Supercodes) Chain Algorithm II, we have 
checked that each of RCi and RC2 is a subcode of any of the five Type II [32, 16, 8] codes. We 
denote the five codes based on RCi {RC2, respectively) by C8l\ . . . , C85^ {C8l'^, . . . , C85'^, 
respectively). 

Hence we obtain: 

Theorem 5.5. Each code C of the five Type II [32, 16, 8] codes has 

ODP[C] = [8, 8, 8, 8, 8, 12, 12, 12, 12, 12, 16, 16, 16, 16, 16, 32]. 

One generator matrix for each Type II [32, 16, 8] code with respect to the ODP in the 
dictionary order is given in the appendix. 



RCi 



11111111111111111111111111111111 



00000000000000001111111111111111 
00000000111111110000000011111111 
00001111000011110000111100001111 
00110011001100110011001100110011 
01010101010101010101010101010101 



10000001000101110100110100100100 
01000001000101000010011110001101 
00100001010001110111010000010010 
00001001000010010101110010100011 
00100001000100100001110111010001 



RC2 



11111111111111111111111111111111 



00000000000000001111111111111111 
00000000111111110000000011111111 
00001111000011110000111100001111 
00110011001100110011001100110011 
01010101010101010101010101010101 



10000001000101110100110100100100 
01000001000101000010011110001101 
00100001010001110111010000010010 
00001001000010010101110010100011 
00100001000100100111101101001000 



Corollary 5.6. For each extremal Type II length 32 code, there are maximum dimension 
subcodes with respect to d = 12, 16, 32 that are involved in the subcode chain for the ODP 
in dictionary order or the inverse order. Furthermore, each extremal Type II length 32 
code contains dimension optimal (and minimum distance optimal) subcodes with parameters 

[32,11, 12], [32, 6, 16], [32, 1,32]. 

5.5 n = 48 

The extended QR code 543 is a unique [48, 24, 12] self-dual code. Using Random (Subcodes) 
Algorithm I, we find that for d' = 16, there is a maximal [48, 14, 16] subcode of q^g. The 
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best known minimum distance optimal [48, 14] code has d = 16. (Note that 17 is the upper 
bound.) One code is given in Magma. We have checked that our code is not equivalent 
to this code. Similarly, for d' = 20, there is a maximal [48,9,20] subcode of q^s- This is 
minimum distance optimal. One [48, 9, 20] code is given in Magma. We have checked that 
our [48, 9, 20] code is not equivalent to this code. For d' = 24, there is a maximal [48, 6, 24] 
subcode of ^43, which is in fact a unique code by [27]. This is minimum distance optimal. 
One code is given in Magma. We have checked that our code is equivalent to this code. 

With respect to the inverse dictionary order we have examined some self-complementary 
subcodes of ^48. There is a [48,5,24] self-complementary subcode (note that A; = 5 is the 
maximum dimension of a [48, k, 24] self-complementary subcode since the unique [48, 6, 24] 
code does not contain the all-one vector). There is a maximal [48, 9, 20] self-complementary 
subcode containing the [48, 5, 24] code (note that /c = 10 is the maximum dimension of a 
[48, /c, 20] self-complementary subcode). 

Lemma 5.7. ([TTl the MacWilliams Identities, p. 129]) Let C he an [n, fc] code and denote 
A^ and A^ to he the numher of codewords of weight w in the code C and C^ respectively. 
Then 



E 



AiPJn, i) = 2''^i, for I) < w < n, 







where Pui{n,i) = J27=oi~^y ( ■ ) ( . \ is a Krawtchouk polynomial. 

Let C be an [n, fc, d] code over Fg. Let T be a set of t coordinates. Let C{T) be the set 
of codewords of C which are on T. We puncture C{T) on T to get a linear code of length 
n — t called the code shortened on T and denoted by Ct |12j . 

Lemma 5.8. ( |12[ Theorem 1.5.7 (iii)]) Let C he an [n. A;, d] code over F^. Let T he a set of 
t coordinates. If t = d and T is the set of coordinates where a minimum weight codeword is 
non-zero, then {C-^)t has dimension n — d — k + 1. 

Both Lemma 15.71 and Lemma 15.81 are useful in determining the non-existence of codes 
with particular parameters and restricted weight distributions. These lemmas are invoked 
to prove the non-existence of particular subcodes of the extended quadratic residue code: 
q^s- Lemma ET] is also applied to determine the possible weight distribution of a putative 
subcode. 

In what follows, we classify all possible weight distributions of a supposed [48, 10, 20] 
self- complementary subcode of ^43. 

Lemma 5.9. If C is a self- complementary [48,10,20] suhcode of q^s, then the non-zero 
codewords of C have weights 20,24,28,48. 

Proof. Suppose to the contrary that C has non-zero weights 20,28,48. Then clearly A20 '■ = 
2^ — 1. Using the MacWilliams Identities (Lemma 15. 7p we obtain the equation 2256+16^420 = 
2^°A^. Hence A^ = ^, a contradiction. D 

Lemma 5.10. If C is a self- complementary [48,10,20] suhcode of q^s, then d-^{C) 7^ 2. 
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Proof. Suppose to the contrary that d^{C) = 2. Shortening C on a minimum weight 
codeword X2 of C^ yields a [46,9,20] code C^q with possible non-zero weights 20,24,28 by 
Lemma [5.81 (here we switched the role of C and C^). 
Define the following matrices: 



B 

A 



K(C46) 



^2o(C*46) 



A24{C,e) A2s{C,,)f. 



Then the MacWilliams Identities yield the matrix equation 2^B = PA, where 



1 

46 

1035 

15180 



1 1 

6 -2 

-5 -21 

-100 44 



1 

-10 

27 

60 



By Grassl's table [9] there (respectively) does not exist a [45,9,20] linear code and there 
does not exist a [44,8,20] linear code, therefore respectively we have AJ^^C^q) = and 
^^(6*46) = 0. Combined with the fact that Aq{C4q) = 1 the above matrix equation yields a 
unique solution of: 

A = [1 243 147 121]^. (1) 

The possible weight distribution of C^q and C^ follows from ([T]). In particular, d{C^ 



46) 



which by shortening C46 on a minimum weight codeword of C^ using Lemma 15.81 implies 
the existence of a [43, 7, 20] code with non-zero weights 20,24,28. This is a contradiction to 
the classification of [43,7,20] due to Bouyuklieva and Jaffe [3j. D 

Lemma 5.11. If C is a self- complementary [48,10,20] subcode of q^s, then there is one 
possible weight distribution of C : 



Ao = l ^20 = 348 A24 
Proof. Define the following matrices: 



326 A28 = 348 ^48 = 1. 



B 
A 



[A^ 
[Ao 



Ai 



AiV, 



A20 A 



24 



Then the MacWilliams Identities along with the fact that C is self-complementary yield the 
matrix equation 2^^B = PA, where 



P 



2 

2256 
389160 



2 1 

16 -24 

-600 276 



By the previous lemma A^ = 0, combined with the fact that Aq = Aq = 1 the above matrix 
equation yields a unique solution of: 

A 



[1 348 326]^. 



D 
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Lemma 5.12. There does not exist a self- complementary [48,k,16] subcode C of q^s for 
k> 17. 

Proof. Suppose a [48,17,16] self-coinpleinentary subcode C exists. The possible non-zero 
weights of C are 16,20,24,28,32,48. Define the following matrices: 



B 

A 






Ai 



Ai]^^ 



Aie A20 A2aY. 



Then the MacWilliams Identities along with the fact that C is self-complementary yield the 
matrix equation 2^^S = PA, where 



2 


2 


2 


1 


2256 


208 


16 


-24 


389160 


40 


-600 


276 


24543024 


-14544 


5616 


-2024 



Isolating the matrix A yields the matrix equation 2^^ P ^B = A where 



21^P- 



17/14 65/224 3/56 1/224 

9729/2 17457/32 211/8 -15/32 

207552/7 -1012/7 -752/7 12/7 

62040 -1605/2 162 -5/2 



The first row of 2^^P ^ implies 

65^ 3^ '^ .^ _ 3 
224 ^ 56 ^ 224 ^ 14' 

which is impossible as Aj- > for all i. Hence no such code C can exist. 



D 



The previous lemmas and example from this section yield the following theorem towards 
the inverse dictionary order ODP for q^g. 

Theorem 5.13. 

OL'P""[g48] = [12, 12, 12, 12, 12, 12, 12, 12, ai, 02, 03, 04, 05, ag, b, 20, 20, 20, 20, 24, 24, 24, 24, 48] 
where ai G {12, 16} and b G {12, 16, 20}. 

Proof. Since ^48 contains the all-one vector, the repetition code [48, 1,48] must be the one 
dimensional subcode first appearing in the subcode chain. By [27] there is a unique [48, 6, 24] 
code with non-zero weights 24, 32; since this code does not contain the all-one vector it cannot 
be involved in the inverse dictionary order subcode chain. Hence k < 5 for a [48, k, 24] 
code involved in the subcode chain. Applying Random (Supercode) Algorithm II to the 
[48, 1, 48] subcode of ^43 we obtained a subcode chain involving a [48, 5, 24] code contained in 
a [48, 9, 20] subcode of g48- Therefore ODF'''"[q^s\i = 24 for 2 < z < 5, and ODP^™[g48]j = 20 
for 6 < i < 9. The maximum dimension for a [48, k, 20] code is fc = 10 by Grassl's table [9], 
hence ODP*™[g48]io = & for 6 G {12, 16, 20} and also ODP^™[g48]j = ^i for 11 < j < 16 and 
tti G {12, 16}. Finally, ODP*""[g48]i = 12 for 17 < i < 24 by Lemma [512 □ 
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Lemma 5.14. There does not exist a [48,k,16] subcode C of q^s for k > 17. 

Proof. Suppose a [48,17,16] subcode of C exists. Since the self-complementary case is already 
considered in Lemma 15.121 we only need to examine the case where the maximum weight 
in C is 36 since the non-zero weights in q^g are 12, 16, 20, 24, 28, 32, 36, 48. Hence the 
possible non-zero weights of C are 16,20,24,28,32,36. Define the following matrices: 



B 
A 



[Aq A;^ A^ A3 A^ A] 

[Ao A16 A20 A24 A28 A32 A 



Air, 



IQ ^16 ^20 ^24 ^28 ^32 ^36j 

Then the MacWilliams Identities yield the matrix equation 2^^B = PA, where 



P 



1 

48 

1128 

17296 

194580 



1 

16 
104 
304 

20 



1 



1712304 -2672 
12271512 -7272 



-104 
-300 
456 
2808 -2024 



1 


-24 


276 




Isolating the matrix A yields the matrix equation 2 P B = A where 



104 
-300 
-456 
2808 

17P-1, 



1 

-16 

104 

-304 

20 



1 

-24 

264 

-1736 

7380 



2672 -19800 
-7272 25080 



2l7p- 



The first row of 2 P implies 



34/21 17/21 65/168 1/6 

4788 1698 2109/4 135 

30000 4592 -61 -312 

61360 680 -965 140 

212448/7 -39488/7 158/7 96 

4482 -1239 3633/8 -81/2 

272/3 -272/3 65/3 -56/3 

17p-l 



1/14 

23 

-92 

132 

-536/7 

19/2 

4 



17 



l-L 



65 



—At + At + -At + —At + —At + At 

21 ^168 ^6 ^14 ^42 ^168 ^ 



1/42 1/168 

1 -3/4 

-8 3 

20 -5 

-160/7 30/7 

25/2 -15/8 

-8/3 1/3 



13 

'21' 



which is impossible as Aj- > for all i. Hence no such code C can exist. 
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Theorem 5.15. 

ODP'^'^iq^s] = [12, 12, 12, 12, 12, 12, 12, 12, a^, a^, 16, 16, 16, 16, h, 62, h, h, Ci, C2, C3, C4, d, e] 

where a^ G {12,16}, bk G {16,20}, q G {16,20,24}, d G {16,20,24,28,32}, and e G 
{20,24,28,32,36,48}. 

Proof. Note that the non-zero weights in ^43 are 12,16,20,24,28,32,36,48. Therefore by 
Grassl's table [9] we may deduce the following: 



ODF [q4s]j = 12 for 17 < j < 24, by the previous lemma. 



ODP'^^"[g48]i = ai for 15 < j < 16 and a^ G {12, 16}. 

ODP *'^[q'48]j = 16 for 11 < j < 14, because as mentioned at the beginning of this sec- 
tion, there exists a maximal [48, 14, 16] subcode of ^48- 

ODP'^^'=[g48]j- = hi for 7 < j < 10 and h E {16, 20}. 

ODP'^^'=[g48]j = Q for 3 < J < 6 and q G {16, 20, 24}. 

ODP'^^"[g48]2 = rf for rf G {16, 20, 24, 28, 32}. 

ODP'^^'=[g48]i = e for e G {20, 24, 28, 32, 36, 48}. 

Note that 16 is not present for values of e because if so then the [48, 14, 16] code involved 
in the subcode chain would have to be a constant weight code. There does not exist a con- 
stant weight code (with weight 16) of dimension greater than 5 by the following reasoning. 

Suppose there exists a [48, k, 16] constant weight code. Define the following matrices: 

B = [A^ Air, 

A = [Ao A,eV. 
Then the MacWilliams Identities yield the matrix equation 2^B = PA, where 

P- 

Since Aq = 1 = Aq, then the matrix equation yields the following system: 



1 1 
48 16 



2^= = 1 + Ai6 
2M^ = 48 + 16^16. 

Solving for Aiq in the first equation and substituting into the second equation yields: 

2M^ = 48 + 16(2*^-1). 
Solving for Aj^ we obtain: 

Ai = 2^-^ + 16. 
And finally since A^ is an integer, then k < 5. D 

From the previous ODPs that have been found for Type II codes, dimension optimal 
subcodes are involved in subcode chains. Therefore we have the following: 

Conjecture: A [48, 6, 24] code is involved in a subcode chain for the ODP in dictionary 
order. 

Thus we have the following corollary. 
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Corollary 5.16. If a [48, 6, 24] code is involved in a subcode chain for the ODP in dictionary 
order, then 

ODP^'^lq^s] = [12, 12, 12, 12, 12, 12, 12, 12, ai, aa, 16, 16, 16, 16, h, 62, h, h, 24, 24, 24, 24, 32, 32] 

where a^ G {12, 16} and bj G {16, 20}. 

We were able to find a doubly-even self-complementary [48, 16, 16] code with generator 
matrix ^[48^16,16]- Such a code was previously not known to exist. Only one singly-even 
self-complementary [48, 16, 16] code was found by A. Kohnert [T5] . 

The dual code has minimum distance d = 4. The generator matrixfor this doubly-even 
self-complementary [48, 16, 16] code is the following: 

100100000000001000110001001011100011100100010100 
010100000000000001111000001001001100110010110100 
001100000000001000101111000010101001010111011110 
000010000000001000101000111100101110101111000111 
000001000000000000011110010000101010101010110101 
000000100000001001001110101100000101100100001101 
000000010000000001111111110010100000011110000000 
000000001000000001100001101101000111111110000000 
000000000100000000110011100111100100110111111100 
000000000010000001001101100111100011001111111010 
000000000001001001110011111110000010100000101000 
000000000000101000010100011110000011000111001110 
000000000000011001111001100000000110000110011110 
000000000000000100001101011001101101011010011000 
000000000000000010010010111001101101000011100110 
000000000000000000000000000000011111111111111110 



G 



[48,16,16] 



Open Problem 1: Determine if the code with generator matrix G'[48,i6,i6] is equivalent 
to a subcode of q^s- 

5.6 n = 72 

Note that qr2 (the extended quadratic residue code of length 72) is a Type II [72, 36, 12] 
code. Due to the complexity, we use Random (Subcodes) Algorithm I. For d' = 16, there is 
a maximal [72, 29, 16] subcode of ^72 with Aiq = 2160. The best known minimum distance 
optimal [72, 29] code has d = 16 (and at most d < 21) with Aiq = 28417, given in Magma. 
Hence our code is not equivalent to this code. For d' = 20, there is a maximal [72, 23, 20] 
subcode with A20 = 3046. The best known minimum distance optimal [72, 23] code has 
d = 20 (and at most < 24) with A20 = 7120 given in Magma. Hence our code is not 
equivalent to this code. 

We start from a best known linear [72, 31, 20] code, given in Magma. Let Ci be this code 
and let d' = 16 < d = 20. Using Random (Supercode) Algorithm II, we have constructed 
in a few seconds a doubly-even self-orthogonal [72, 35, 16] code C containing Ci with Aiq = 
129972. It is known from Magma that there is a best known minimum distance code with 
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parameters [72, 35, 16]. This is a doubly-even self-orthogonal code with Ai^ = 136116. Hence 
our code is not equivalent to the known code. We do not know how many doubly-even self- 
orthogonal [72, 35, 16] codes exist. 

6 Conclusion 

The optimum distance profile for a linear code (and any code in general) is a relatively new 
concept developed in [3] and [12]. This area is particularly interesting due to its practical 
applications. In this paper we relate the optimum distance profile of a code to the concept of 
maximal subcodes of high minimum distance. We develop four algorithms which are highly 
efficient in comparison to a brute force examination of all subcodes. 

The classification of self-dual codes continues to be an extremely active area in coding 
theory. A particularly interesting class of self-dual codes is those of Type II which have high 
minimum distance (called extremal or near-extremal). It is notable that this class of codes 
contains famous unique codes: the extended Hamming [8, 4, 4] code, the extended Golay 
[24, 12, 8] code, and the extended quadratic residue [48, 24, 12] code. A long standing open 
problem in coding theory is to prove the existence or non-existence of a Type II [72, 36, 16] 
code. The aim of this paper is to shed light on the structure of this interesting class of 
codes. We examine the maximal subcodes and ODPs of Type II codes for lengths up to 32. 
Of recent significance is the classification of length 40 Type II codes [1]. The examination of 
these codes would be extensive work as there are 16470 Type II [40, 20, 8] codes (the highest 
minimum distance in this case is 8 which is not minimum distance optimal by [9j). Therefore 
we examined a more interesting case, the unique Type II code 943 of length 48. 



Appendix 



C8V 



RC 



10000000000000010001011000001110 
01000000000000100001010100110001 
00100000000000100011000101001001 
00010000000000010001000011001101 
00001000000100000010010100110010 



C82i 



RCi 



10000000000000010110111000011111 
01000000000000010010000000111011 
00100000000000010110110111010101 
00010000000000010000111001000101 
00001000000000010001101000010011 



C831 



RCi 



10000001000100010011010101100110 
01000001000100100011011001011001 
00100001000100100001001000100001 
00010001000100010010001000100010 
00001001000000000001011100100010 



CM^ 



RCi 



10000000000000010001011000001110 
01000000000000100001010100110001 
00100000000000100011000101001001 
00010000000000010000000110110101 
00001000000100000011010001001010 
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C85^ 



i?Ci 

10000000000000010001011000001110 
01000000000000010101011111011010 
00100000000000010001101000110100 
00010000000000010000000110110101 
00001000000000010000010010011011 



CSV 



RC2 

10000000000100000010010100101100 
01000000000100110010011011101100 
00100000000100110000001010010100 
00010000000100000010001111101111 
00001000000000010001011011101111 



C82^ 



RCo 



10000001000101110001011101111110 
01000001000101000001010001000001 
00100001000100100001001000100001 
00010001000100010001000100010001 
00001001000001100000011000001001 



C83^ 



RCo 



10000001000100010011010101100110 
01000001000100100011011001011001 
00100001000100100001001000100001 
00010001000100010010001000100010 
00001001000000000001011100100010 



C84^ 



RCo 



10000000000100000010010100101100 
01000000000100110010011011101100 
00100000000100110000001010010100 
00010000000100000011001001101000 
00001000000000010000011101101000 



C85' 



RCo 



10000000000100000010010100101100 
01000000000100000110011111110100 
00100000000100000100001110001100 
00010000000100000010001111101111 
00001000000000010001011011101111 
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